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ABSTRACT 

The response of the ionospheric E-region to solar-geomagnetic storms can be 
characterized using observations of infrared 4.3 um emission. In particular, we utilize 
nighttime TIMED/SABER measurements of broadband 4.3 um limb emission and derive 
a new data product, the NO + (v) volume emission rate, which is our primary observation- 
based quantity for developing an empirical storm-time correction the IRI E-region 
electron density. In this paper we describe our E-region proxy and outline our strategy for 
developing the empirical storm model. In our initial studies, we analyzed a six day storm 
period during the Halloween 2003 event. The results of this analysis are promising and 
suggest that the ap-index is a viable candidate to use as a magnetic driver for our model. 

1. INTRODUCTION 

Observations of thermospheric infrared emission from the TIMED/SABER instrument 
have fostered the development of new data products, models, and analysis tools for the 
study of upper atmospheric and ionospheric response to solar-geomagnetic disturbances. 
The focus of this paper is on analyzing SABER 4.3 um limb emission measurements to 
quantify the E-region response to solar-geomagnetic storms, with the objective of 
developing an empirical E-region storm-time correction to the International Reference 
Ionosphere (IRI) model (Bilitza, 2001). 

The Sounding of the Atmosphere using Broadband Emission Radiometry (SABER) 
instrument is a broadband infrared limb sounder on the Thermosphere-Ionosphere- 
Mesosphere-Energetics and Dynamics (TIMED) satellite (Russell et al., 1999). TIMED 
was launched in 2001 and has observed the atmospheric response to several major solar 
eruptive events. Specific to this paper, SABER observed enhancements in the nighttime 
4.3 um limb emission channel by several orders of magnitude during the April 2002 and 
October-November (“Halloween”) 2003 storm period, for example. The storm-time 
enhancements in the nighttime 4.3 um radiances are due to emission from vibrationally 
excited NO + , i.e. NO + (v) (Mertens et al., 2007). Precipitating auroral electrons during 
solar-geomagnetic disturbances increase the production of ions via ionization of the 
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neutral atmosphere (Banks et al., 1974; Strickland et al., 1976). In the E-region, these 
ions react with neutral species to produce NO + (Torr et al., 1990; Fox and Sung, 2001). 
Some of the ion-neutral reactions are exothermic enough to produce NO + (v), which emits 
at 4.3 um (Winick et al., 1987). The NO + (v) component of the total 4.3 um emission is 
extracted from the SABER measurements by subtracting off the CO 2 CV 3 ) contribution 
using non-LTE radiation transfer algorithms and input atmospheric data from SABER 
measurements, combined with the NRLMSIS-00 model (Mertens et al., 2007, 2002, 
2001). Since NO + is the terminal E-region ion, NO + (v) 4.3 um emission is an excellent 
proxy to characterize the response of the E-region to solar-geomagnetic storms. 

Moreover, by charge neutrality, NO + (v) emission is also an indication of the response of 
the E-region electron density to magnetic storms. We will exploit these properties of the 
E-region in our development of an empirical storm-time correction to the IRI model E- 
region electron density. 

The IRI model is a widely used empirical model for the specification of ionospheric 
parameters and is recommended for international use by the Committee on Space 
Research (COSPAR) and the International Union of Radio Science (URSI). However, the 
specification of the ionospheric response to solar-geomagnetic disturbances in IRI 
remains largely incomplete, and there is currently no storm-time correction to IRI 
parameters in the E-region. An error of several orders of magnitude in the IRI E-region 
parameterization will limit the model’s usability in radio wave propagation models - as 
applied to system design of communication, navigation, and surveillance systems - and 
real-time processing of the radio wave propagation data during solar-geomagnetic storms. 
Errors in E-region electron density of several orders of magnitude can induce non- 
negligible errors in modeling vertical and slant path radio wave propagation (satellite to 
ground communication). Moreover, significant errors can be induced in modeling 
satellite to satellite radio wave propagation along long horizontal paths. To improve the 
utility and range of applicability of the IRI model, we have initiated a plan to develop an 
empirical storm-time correction to the IRI E-region electron density using SABER 4.3 
um limb radiance measurements. 

In this paper we describe the parameterization of the storm-time correction to the IRI E- 
region electron density using observations from TIMED/SABER 4.3 um limb emission 
measurements. Our E-region proxy derived from the SABER observations is the NO + (v) 
volume emission rate (VER), which is described in section 2. The parameterization of E- 
region electron density enhancements due to solar-geomagnetic forcing, based on the 
NO + (v) VER, is described in section 3. Analysis of several magnetically disturbed days 
during the Halloween 2003 storm period, as presented in section 4, provide considerable 
guidance on the selection of ( 1 ) the external magnetic driver index , and ( 2 ) the form of 
the impulse-response function, both of which are key ingredients in the formulation of the 
parameterization of the storm-time correction to the IRI E-region electron density. A 
summary of our initial studies is presented in section 5. 

2. E-REGION PROXY: NO + (v) VER 
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As discussed in the introduction, NO + (v) 4.3 um emission is an excellent proxy for 
characterizing the response of the E-region to solar-geomagnetic disturbances. 
Specifically, we derive NO + (v) VER from the SABER 4.3 um limb emission 
measurements. 

The first step to derive NO + (v) VER is to subtract the background CO 2 CV 3 ) emission from 
the measured 4.3 um limb radiance. This is accomplished with the radiation transfer 
algorithms used in the operational processing of the SABER data. The CO 2 CV 3 ) 
contribution is removed using the 4.3 um forward radiance component of the operational 
non-local thermodynamic equilibrium (non-LTE) kinetic temperature (Tk) and CO 2 
volume mixing ratio (vmr) retrieval algorithm (Mertens et al., 2001, 2002). The forward 
model is comprised of two parts: (1) the CO 2 vibrational temperature (Tv) model and (2) 
the limb radiance model. Limb radiance is calculated using BANDPAK (Marshall et al., 
1994), which is based on emissivity databases calculated line-by-line using LINEPAK 
routines (Gordley et al., 1994). The CO 2 Tv model is based on the Modified Curtis 
Matrix approach (Lopez-Puertas et al., 1986a-b, 1998), and uses BANDPAK in all the 
radiation transfer calculations. 

The nighttime CO 2 CV 3 ) contribution to the SABER-measured 4.3 um limb emission is 
simulated with the non-LTE radiation transfer algorithms described above using input 
atmospheric data from SABER-retrieved Tk and pressure, CO 2 vmr from a TIME-GCM 
climatological database (Roble et al., 1995), and composition data (N 2 , O 2 , O) from the 
NRLMSIS-00 model (Picone et al., 2002). All input composition data (CO 2 , N 2 , O 2 , O) to 
the CO 2 Tv model for the nighttime 4.3 um limb radiance calculations are consistent with 
the operational processing of SABER data below 120 km. Above 120 km Tk and 
composition are obtained from the NRLMSIS-00 model (CO 2 still comes from the TIME- 
GCM database). Pressure is extended above 120 km assuming hydrostatic equilibrium. 

The OH(v) VER data product from the SABER 2.0 um channel measurement is used to 
calculate the vibrational excitation of N 2 U) from collisions of ground state N 2 with 
OH(v < 9). The method of calculating the total production rate of OH(v) from measured 
OH(v) VER at 2.0 um is by Mlynczak et al. (1998), while the method of calculating the 
rate of production of N 2 (l) from OH(v) is by Lopez-Puertas et al. (2004). Energy transfer 
from vibrationally excited OH (i.e., OH(v < 9)) to N 2 (l), and then to CO 2 CV 3 ) via a V-V 
transfer mechanism, is an important source of nighttime CO 2 4.3 um emission starting 
from 70 km and extending to above 120 km (Lopez-Puertas et al., 2004). The influence 
of this excitation process extends well above the OH layer due to the non-LTE process of 
local excitation (OH(v<9) -> N 2 O) -> CO 2 CV 3 )) followed by radiative loss and re- 
absorption at higher altitudes. 

Once the nighttime CO 2 CV 3 ) 4.3 um limb radiance is calculated, it is subtracted from the 
SABER 4.3 um radiance measurement. The residual radiance is emission from NO + (v). 
The second step to deriving the NO + (v) VER is to perform an Abel inversion on the 
residual radiance. An example of the NO + (v) VER from the April 2002 and Halloween 
2003 storm periods is shown in Figure 1. The significant enhancement of NO + (v) VER 
during the storm periods is clearly evident. Note also that the peak in the storm-enhanced 
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NO + (v) occurs roughly between ~ 100-110 km, which is consistent with the altitude of 
the peak in the E-region electron density. 


SABER— Derived N0 + (v) 4,3 um Volume Emission Rate (VER) 


180 

170 


160 


150 

? 

4- 140 

( L 

1 130 
< 

120 


110 


100 
90 

-0.2 0.0 0.2 0.4 0.6 0.8 1.0 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 

VER (ergs cm -3 s' 1 ) x 10' VER (ergs cm" 1 s' 1 ) x 10' 


April 2005 


October 2003 




Figure 1: SABER-derived NO + (v) VER for representative days prior to and during the April 2002 
and Halloween 2003 solar storms. The blue lines correspond to NO + (v) VER prior to the storms, 
while the red lines correspond to NO + (v) VER during the storms. Figure taken from Mertens et al. 
(2007). 


3. EMPIRICAL MODEL PARAMETERIZATION 

The empirical storm-time correction to the IRI E-region electron density is based upon 
deriving a storm-induced correction factor from the SABER-derived NO + (v) VER, and 
developing a fit of the correction factor to an appropriate external magnetic driver index. 
The motivation for this approach is that it has been successfully applied to characterizing 
the response of a number of geospace parameters to external drivers (e.g., Fuller-Rowell 
et al., 2000; Vassiliadis et al., 2002). 

The correction factor is defined as the ratio of the storm-time NO + (v) VER to a quiet- 
time, climatological NO + (v) VER, such that the correction factor is given by (Mertens et 
al., 2007) 


r(z.K.<P„.t)=[VER(z,K„<P„.t) ] s l[VER(z,l m ,<p m )\ 


J Quiet 
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In the above equation, VER refers to NO + (v) VER, X m is magnetic latitude, and cp m is 
magnetic local time. The quiet-time, climatological NO + (v) VER are monthly-averaged 
quantities. Results from SABER 4.3 um radiance measurements, NOAA/POES auroral 
electron energy characteristics in-situ observations, charge neutrality, and model 
simulations of the ionospheric plasma, NO + (v) chemical and radiative energetics, 
suggests the NO + (v) VER ratio in (1) is approximated by 
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where [e] is the electron density. The NO + (v) VER ratio will be computed for all 
magnetically disturbed periods contained in the SABER dataset from 2002-2007. 

The computed storm-time correction factors in (1) are fit to a polynomial expansion in a 
magnetic forcing parameter (G), such that 

N 

y 'j (A \n i Rm (.Z, A m , (p m , t ), (3) 

i=0 

where a t are the expansion coefficients. The magnetic forcing parameter (G) captures the 
dynamical response of the E-region by representing the forcing parameter as a 
convolution of an external magnetic driver index (D) with a dynamical response function 
( F(t) ). Thus, the forcing parameter is given by 

T 

G(z , A m ,<p m -,t)= J F{z , Xn ,ip m ]r)D(t — r)dr. (4) 

-T, 

The effective memory of the storm-time response of the E-region corresponds to time (T) 
in the upper limit of the integral in (4). The lower limit of the integral (Ts) represents the 
start time and is chosen to ensure numerical stability in solving for the response function 
(F(t) ). Both upper and lower integration limits can be established by examining the 
cross correlation function presented in section 4. 

To develop the empirical storm-time correction model, one must select an appropriate 
external magnetic driver index (D). To be of practical value to the broad space weather 
community, and for easily implementation into the IRI model, the external magnetic 
driver index must be a widely known, easily accessible index. Some obvious indices to 
choose from are the Hemispheric Power index (HP), Disturbed Storm Time (Dst) index, 
ap index, and the Auroral Electrojet (AE) index. The selection of the appropriate 
magnetic driver index is discussed further in the next section. On the other hand, once the 
magnetic driver index is chosen, the response function can be readily determined by 
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r T 

r(z,\n,(p m ,t) = I F{z,\n,ip m \T)D{t-T)dT. 
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The above equation assumes a linear impulse-response relation between the external 
driver index and the response of the E-region ionosphere, as characterized by the NO + (v) 
VER ratio. The NO + (v) VER ratios for all magnetically disturbed periods in the SABER 
dataset from 2002-2007 are used to solve for the response function (F(r)) in (5). The 
integral in (5) is discretized and formulated in terms of a matrix-vector equation. The 
response function is determined by inverting the equivalent matrix-vector equation using 
Singular Value Decomposition (SVD) techniques (Press et al., 1992). Any residual non- 
linearity between the driver (D) and the response (r) is accounted for by non-linear terms 
in the polynomial expansion in (3). 

4. ANALYSIS OF THE HALLOWEEN 2003 STORM 

In this section we analyze six days during the Halloween storm period, October 27- 
November 1, 2003. The data volume based on six days is insufficient to develop the 
statistical parameterization presented in the previous section, yet this limited analysis 
provides considerable guidance into the selection of magnetic driver and insight into the 
form of the response function. 

To obtain a global perspective on the response of the E-region during the Halloween 
2003 storm period we show the SABER-derived NO + (v) VER at 1 10 km in Figure 2. 
Large enhancements are noticeable on all of the days shown in Figure 2. The morphology 
of the peak NO + ( v) VER is consistent with the shape of the auroral oval due to electron 
precipitation (not shown). Exceptional equatorward expansions of the auroral oval are 
clearly evident on 29-30 October. 

Guidance into the selection of the magnetic driver index can be obtained by computing 
the cross correlation function between the SABER-derived NO + (v) VER and the 
candidate magnetic driver index. The cross correlation function is defined as (Vassiliadis, 
2002) 

T 

G(A m ,r) = 1 i f2( VER (\n,t t + r) - ( VER(X m )))(D(t t ) - (D)\ (6) 

1 cr D a V ER v x ' ' x ' 

The data volume for the six day analysis period (October 27-November 1, 2003) is 
insufficient to determine a dependence on magnetic local time. Thus, the NO + (v) VER 
data during the six day period were sorted and averaged into 5 degree bins in magnetic 
latitude. Both the NO + (v) VER and the magnetic driver indices were sorted an averaged 
into 3 hour UT-time bins throughout the duration of the six day storm period. The 
summation in (6) is over all the 3 hour UT-time bins over the six day storm period. The 
average values (denoted by <>) were computed over the six day period and the standard 
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deviation (denoted by a ) were computed about the average values. 



Nighttime NO + (v) Volume Emission Rote, Halloween Storm 2003 at Alt=110km 
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Figure 2: : SABER-derived NO + (v) YER at 110 km during a select period of the Halloween 2003 
solar storm event. The unit for the VER is ergs cm' 3 s' 1 . 

Figure 3 shows the cross correlation function between NO + (v) VER and the Hemispheric 
Power (denoted HP: unit is GW) index derived from NOAA/POES measurements of 
auroral electron precipitation. Since the HP-index most likely represents to the physical 
source mechanism for the response of the E-region to solar-geomagnetic storms, we 
consider this correlation to be a benchmark for evaluating the “goodness” of more readily 
available solar-geomagnetic storm indices. Figure 3 indicates that the peak correlation of 
~ 0.60 for altitudes at and below 1 10 km occurs between 50 and 60 degrees magnetic 
latitude, and occurs instantaneously or with delay times of the order of ~ 3 hours UT. At 
and above 110 km the peak value of the cross correlation extends down to ~ 30-35 
degrees magnetic latitude, and tilts to the right with increasing delay time. This is likely 
an indication of composition changes in the region of auroral electron dosing, which are 
then transported equatorward. 

We have computed the cross correlation function between NO + (v) VER and the Dst- 
index (not shown). The Dst-index is indicative of perturbations of the horizontal 
component of the geomagnetic field at the magnetic equator due to storm-time 
disturbances of the ring current. The perturbation in the horizontal magnetic field 
component could potentially alter the mapping of the auroral particles into the 
atmosphere at mid- and high-latitude, and subsequently influence the response of the E- 
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region to auroral particle precipitation. The cross correlation between NO + (v) VER and 
Dst-index has a similar morphology as shown in Figure 3, but the magnitude of the 
correlation as far less than between NO + (v) VER and HP-index. 
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Figure 3: Cross correlation between SABER-derived NO + (v) VER and Hemispheric Power (denoted 
HP: unit is GW) index derived from NO A A/POES satellites. Each panel shows the cross correlation 
function as a function of magnetic latitude and delay time (hours). The different panels correspond to 
different altitudes. The cross correlation function is defined in the text. 

Figure 4 shows the cross correlation between NO + (v) VER and the ap-index. The 
correlation between these two quantities is comparable to the correlation between NO + (v) 
VER and the HP-index. This is a fortunate result because the ap-index is readily 
accessible and the IRI model already includes the ap-index as an input. We haven’t 
examined the AE-index yet, but we expect this index to be very useful as well, 
particularly at the higher latitudes. 

Figure 5 shows the solution for the response function based on the six days analyzed 
during the Halloween storm period. The results are very noise due to the small data 
volume. However, a reasonable pattern can be discerned which is consistent with the 
cross correlation functions shown in Figures 3 and 4. Poleward of 50 degrees magnetic 
latitude, where the largest enhancements occur, a peak in the response functions occurs 
between zero and a ~ 3 hour FIT delay time. This is consistent with the delay times 
indicated by the cross correlation functions in Figure 3 and 4. For larger delay time the 
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response function oscillates through zero and eventually drifts positive. We expect that 
when we analyze all the magnetic storm periods in the SABER database from 2002-2007 
the peak in the response function at zero to 3 hour UT delay time will dominate and the 
response function will quickly decay to zero with small oscillations about zero. 
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Figure 4: Cross correlation function between SABER-derived NO + (v) VER and the ap-index. Each 
panel shows the cross correlation function as a function of magnetic latitude and delay time (hours). 
The different panels correspond to different altitudes. The cross correlation function is defined in the 
text. 


5. SUMMARY 

The NO + (v) VER derived from SABER 4.3 um limb emission measurements provide an 
excellent proxy for characterizing the response of the ionospheric E-region to solar- 
geomagnetic storms. Storm-time observations from SABER and NOAA/POES, 
combined with the known physics and chemistry of the E-region, suggest that the NO + (v) 
VER also contains valuable information about the response of the E-region electron 
density to solar-geomagnetic disturbances. Thus, we have initiated a program to develop 
an empirical storm- time correction to the IRI E-region electron density based on NO + (v) 
VER derived from SABER. Our storm model will potentially enhance the capability of 
the widely used IRI model and increase its usability in real-time space weather 
applications. Preliminary analysis of a six day storm period during the Halloween 2003 
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event suggests that the ap-index is a viable candidate for an external magnetic driver for 
our empirical model. Additional storm indices will be examined in the future. Moreover, 
we will analyze over 80 storm days observed in the SABER dataset from 2002-2007. The 
results of this analysis will provide the input to develop the final form of our empirical 
storm model, which will be the subject of future reports. 
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Figure 5: Impulse-response function as a function of delay time (hours) derived from a six day period 
during the Halloween 2003 storm event. The figure shows the response function derived at several 
magnetic latitudes. The response function is defined in the text. 
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